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Abstract. In this paper we provide a novel strategy to prove the validity 
of Hartree's theory for the ground state energy of bosonic quantum systems 
in the mean-field regime. For the known case of trapped Bose gases, this 
can be shown using the strong quantum de Finetti theorem, which gives 
the structure of infinite hierarchies of fc-particles density matrices. Here we 
deal with the case where some particles are allowed to escape to infinity, 
leading to a lack of compactness. Our approach is based on two ingredients: 
(1) a weak version of the quantum de Finetti theorem, and (2) geometric 
techniques for many-body systems. Our strategy does not rely on any special 
property of the interaction between the particles. In particular, our results 
cover those of Benguria-Lieb and Lieb-Yau for, respectively, bosonic atoms 
and boson stars. 
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1. Introduction 

In this paper, we consider a quantum system composed of a very large number 
of interacting particles. Because of the correspondingly large number of degrees 
of freedom, it is extremely hard to describe the precise behavior of a system of 
this kind. It is thus often useful to resort to approximate theories which are 
simpler to deal with. One of the major issues in many-body physics is then 
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to justify the validity of the effective models, that is, to relate them to the 
many-body problem in a particular regime. 

For bosons, the simplest effective theory can be obtained by assuming that all 
particles are in the same quantum state. This leads to the celebrated nonlinear 
model introduced by Hartree in |28j Q In this theory, the particles in the system 
behave as if they were independent, but submitted to a common mean-field 
potential due to all the other particles. 

The purpose of this paper is to prove that Hartree's theory gives a correct 
approximation of the ground state of the many-body bosonic system in the limit 
of large particle number TV with the intensity of the pair interaction assumed to 
decrease proportionally to A -1 (which is often called the mean-field regime). 
As we will recall below, there are many results of this kind in the literature, 
most of them dealing with particular systems. Here we propose a novel method 
which allows to deal with a very large class of many-body systems and does 
not depend on the special form of the interactions. We will be particularly 
interested in the case where some particles are allowed to escape to infinity, 
leading to a possible lack of compactness. 

The model. Let us consider a system composed of N identical bosons. The 
one-particle space is any separable Hilbert space S) and the whole system is 
therefore described by the A-fold symmetric tensor product fi N := (££) s $). We 
assume that the TV-body Hamiltonian takes the following form: 

N 1 

Hn:= Y j T 3 + w — i £ w M- (1.1) 

j=l l<k<e<N 

Here T is a self-adjoint operator on S) that accounts for the self energy of 
the particles, and w is a symmetric operator on the two-particle space fj 2 , 
which corresponds to the pair interactions between the particles. As usual, 
Tj = 1®...®T<8>...®1 denotes the associated operator acting on the j-th 
particle and Wkt is the potential acting on the pair (k, £) of particles. We will 
always assume that Hn is bounded from belowjj Note that, even if we stick to 



Hamiltonians of the form (1.1) for simplicity, most of our results will indeed be 
valid for a much larger class of mean-field Hamiltonians. 

The fact that we are considering the mean-field regime is apparent in the 



factor 1/(A — 1) in front of the interaction term in (1.1). It has the effect of 
keeping the single particle energy and the interaction energy of the same order 
of magnitude, so that one may expect a well-defined limit problem. Note that 
this factor could be replaced by any constant behaving like 1/N in the limit 
./V — > oo, without changing the result; the use of 1/(A — 1) only simplifies some 
expressions. While this is certainly not the only scaling one may consider, it 
is simple and instructive, and has been very often considered in the past as 
a model case for the rigorous derivation of mean-field theories in many-body 
physics. 

1 More precisely, the theory was intended to be used for the electrons in an atom but it had 
to be corrected later by Fock since electrons are not bosons, but rather fermions. 

o 

In the mean-field regime, this automatically implies Hn > —CN. 
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For physically relevant examples, we may typically think of bosons living in 
a bounded set Q C M d and T = — A in fj = L 2 (Q) with appropriate boundary 
conditions, or think of particles in M rf and T = —A + V(x) in S) = L 2 (JR. d ). 
In the latter case, either V(x) — > oo as \x\ — > oo (V is then called a trapping 
potential), or V(x) — > as \x\ — > oo but V(x) is negative somewhere to bind 
(some of) the particles. The one-body operator T may also involve the magnetic 
Laplacian (— iV + A(x)) +V(x) corresponding to a given vector potential A on 
M 3 , or a pseudo-relativistic operator \J m 2 c^ — A — mc 2 + V(x). The two-body 
potential w is often the multiplication operator by an even real-valued function 
w(x — y) that decays at infinity. 

The Hartree functional is obtained by restricting the energy functional (qua- 
dratic form) of Hn to uncorrelated functions of the form \& = u® N , where u is 
a normalized vector of Sj. This leads to the nonlinear Hartree energy 



u 



® N H N u® N ) I 



y (u,Tu)% + -(u®u,wu®u) s p :=£n(u). (1.2) 

If E(N) = inf a(Hjf) denotes the bottom of the spectrum (ground state energy) 
of Hn, it is then expected that 

l im ^)=e H (1 .3, 

where en is the minimal Hartree energy 

e H := inf £ u (u). (1.4) 

1 1 u 1 1 = 1 

Since ^ = u® N can be used as a trial state, it is obvious that E(N)N~ 1 < en 
for all N > 2. The lower bound is much more subtle and it means that the 
purely uncorrelated ansatz u® N does capture the first order of the ground state 
energy in the limit iV — > oo. We will come back below to the consequences it 
may have on the ground state of fljy (when it exists), such as Bose-Einstein 
condensation. 

Our aim in this paper is to provide a general strategy to justify the conver- 



gence (1.3), which is applicable to a very large class of physical systems. The 



proof of (1.3) is available in the literature for numerous special cases, including 
"bosonic atoms" [H [62j [5l El [33], boson stars [IS 06], the homogeneous Bose 
gas [56] , trapped Bose gases [27] , the Lieb-Liniger model |57J , and many 
others. More abstract models are discussed in (23] [691 E21 ED] - Recently, the 
experimental realization of Bose-Einstein condensates has motivated the study 
of a more subtle, so called Gross-Pitaevskii, limit |431 141j. We will not consider 
this limit here, but we hope that our method will in the future be useful to deal 
with it as well. 

It is striking that the Hartree approximation is valid in very different physical 
situations. For example, in bosonic atoms the interactions are repulsive and 
the particles are submitted to an attractive potential generated by fixed nuclei, 
whereas in boson stars the interactions are attractive and the system is fully 
translation-invariant. This generality holds despite the fact that most proofs 
rely on specific properties of the Hamiltonian (and in particular of the two-body 
potential w). The main message of the present article is that the validity of the 
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Hartree approximation does actually not rely on any specific properties of the 
Hamiltonian, but is rather a consequence of the special structure of the set of 
bosonic density matrices for large N. 

Before explaining this, let us insist on the fact that we are interested here 
in the large-iV behavior of the ground state energy E(N) of the Hamilton- 
ian. There are many works on the related (but still different) derivation of the 
time-dependent Hartree theory from the time-dependent Schrodinger equation 
associated with the Hamiltonian Hn, see for instance [291 123 [Ml El [201 EH 
[3j [22l [2H [5H [Ml EI] ■ I n this case one starts close to a Hartree state at time 
zero, and then proves that the Schrodinger flow stays close to the correspond- 
ing trajectory of the Hartree state. It is fair to say that the validity of the 
Hartree approximation has been proved under much more general assumptions 
in the time-dependent case than for the ground state energy E(N). Our work 
will therefore place the time-independent problem on the same footing as the 
time-dependent problem. 

Representability and de Finetti theorems. Our starting point is the for- 
mulation of the problem in terms of reduced density matrices [42 . The k- 



(k) AT 

particle density matrix 7^/ of a pure iV-body state 6 Sj is defined by 

7i fe) :=Tr fe+1 ^ Ar |*)(v[>| (1.5) 

for < k < N, where Ti^+i^n denotes the partial trace with respect to the 
last N — k variables and is the orthogonal projection onto the state *. 

(k) h (k) 

Thus 7^, is a positive trace class operator on fj with Tr^k 7^ = 1 (note the 
normalization convention) . 

The main interest of density matrices is that the energy per particle can be 
expressed^] only in terms of 7^ : 

- Tr, (T 7 U) + ilV fa?) =^(H^). (1.6) 



We see that, thanks to the mean-field factor l/(iV— 1) in front of the interaction, 

the expression of the energy is even completely independent of N. The N 

(2) 

dependence is hidden in the constraint that 7^, must arise from an iV-body 
state One may thus reformulate the ground state energy of our system as 

E{N) - Snf^ -^-Prt ~(2)^(2) 



{^{H^^VevP} (1.7) 



N 
where 

= { 7 (2) E © W) :3«€fl w , ||*|| = 1, 7 {2) = 7? } } 

is the set of all the two-particle density matrices arising from a pure iV-body 
state \P ( u N-representable" two-particle density matrices). Here (3 1 (^) is the 
space of all trace-class operators on the Hilbert space 

It is often very useful to work with mixed states instead of pure states. A 
mixed state is an operator G on $j N which is a convex combination of pure 



3 And of course, more generally, a problem involving at most fc-particle interactions only 
depends on the fc-particle density matrix. We stick to the two body case for clarity. 



DERIVATION OF HARTREE'S THEORY 



5 



states: G = Y2i n i\^i)(^i\ with n, > and Yli n i = 1- As the density matri- 
ces 7^ depend linearly on the operator |^ r )(^ r |, their definition can easily be 
extended to mixed states: 

7 g° :=Tr k+1 ^ N G. (1.8) 
The iV-body energy of a mixed state is now 

Tr (H N G) _ ( 1 / (2) \ _ 1 (2)> 



Tr fi (T 7 ^j + 2 TV (w 7 ^j = - ^ a (H 27 ^)- (1.9) 

Diagonalizing G and using that the energy is linear in the two-body density 
matrix, one sees that minimizing over mixed states gives the same answer as if 
one minimizes over pure states only. Therefore, we can also write 

^ = iinf{Tr, 2 (^)), 7 ^ G <} (1.10) 

where 

VP = G ©W) : 30<GG6V), Tr^ G = 1, ^ = ^} 

is the set of mixed-representable two-particle density matrices, which coincides 

~ (2) 

with the convex hull of the set V%' . Of course we may define in a similar fashion 

the sets and of /c-particle density matrices, arising respectively from 
pure and mixed states. 

For fixed N, rewriting the ground state energy as in (1.7) or (1.10) is not 



particularly helpful. On the contrary to the one-particle set V N which is well 

(2) 

understood, it is indeed a famous open problem to characterize the set V N . 
This is called the N -representability problem, usually stated for fermions |13j . 

For bosons, : 
now explain. 

Taking 
sequence 



For bosons, it is possible to describe the sets V N in the limit N — > oo, as we 
Taking partial traces it is easy to see that the sets pjy form a decreasing 



' N+l L ' AT • 



One may then see Problem (1.10) as the minimization of a fixed energy func- 
tional on a variational set that gets more and more constrained as N increases. 
The energy E(N) /N is thus increasing with N and with a leap of faith one may 
hope that our variational problem will converge to the one posed on the limit 

inti 

7>W := p pW (i.n) 



(k) 

of the sequence of sets V N , that is on their intersection: 



N>1 



the set of /c-body density matrices that are iV-representable^] for any N > k. If 
we are allowed to exchange the infimum in ( |1.10[ ) and the limit ( |1.11[ ), we then 



^Remark that the set is empty for fermions since 7^ < 1 and Trj^ = 1 with 



our choice of normalization. 
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formally obtain 



Urn 



inf |Tr^2 



H 2 j 



(2) 



p(2)J 



(1.12) 



Provided we can pass to the limit to write (1.12), the validity of Hartree's 



theory then follows from the fact that is the convex hull of the two-particle 
density matrices of Hartree states, |u® 2 ) (ii® 2 |, as we now explain. 

Describing the structure of the limiting sets is precisely the object of the 
so-called quantum de Finetti theorem, proved by St0rmer and Hudson-Moody 



in |65l I5T] and recalled in Theorem 2.1 below. This result is a quantum gener- 
alization of the famous classical de Finetti, also called Hewitt-Savage, theorem 
|14l [T5l \19\ \30\ [TB] about symmetric probability measures having an infinite 
number of variables. The importance of such results for mean-field theory 
has been known for a long time in the context of classical statistical mechan- 
ics [TD1 El E3 D3J E2]- In [B3 E3 EZ! Varadhan's large deviation principle 
(similar in spirit to the classical de Finetti theorem) was used to understand 
the Bose-Einstein condensation of certain quantum systems. The quantum ver- 
sion of the de Finetti theorem was then used to treat a larger class of problems 
in [231 EQl [521 [70] . 

In our language, the quantum de Finetti theorem simply states that, for any 
fixed k, the limiting set is the convex hull of the k-particle density matrices 
of Hartree states, the latter being its extremal points. The operators G J>^) 
can therefore all be written in the form 



7 



(k) 



,®k\ l„,®k 



){u m \dfi(u) 



where /i is a Borel probability measure on the sphere SSj of the one-particle 



Hilbert space S). Now, we can compute the right side of (1.12): 



1 



Tr.2 tf 27 



,(2) 



IT 



Mm* 



dfj,(u) 



SSj 



£n(u) dfj,(u) > eu 



where in the last inequality it is used that fi is a probability measure. Hence 



the right side of (1.12) is nothing but en- 



Of course, in the argument that we have sketched above, the main difficulty 



is to justify the formal limit (1.12). A typical case that can be easily dealt with 



is that of trapped particles, that is, when the single-particle Hamiltonian T has 
a compact resolvent. Simple examples consist of non-relativistic particles living 
in a bounded domain, or T = — A + V with V a trapping potential. Proceeding 
as sketched above, one can easily justify Hartree's approximation for trapped 
bosons. We quickly study this situation in Section [3j This is very much in 
the spirit of the earlier works |23[ [52] and we include it mainly for pedagogical 
purposes. Interesting cases covered by this approach include the homogeneous 
and trapped Bose gases. 

In many practical cases, however, the particles are not all trapped and some 
can escape to infinity. It should then be verified that those escaping to infinity 
are still correctly described by Hartree's theory. The situation is therefore much 
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more complex and a more detailed analysis is necessary. It is the main object 
of this article to provide a strategy to carry over this detailed analysis. 

In infinite dimensions, one has to be very careful of the topology which is 



used for investigating the limit of V^P- The set defined in (1.11) is actually 



(k) 

the limit of V N for the trace norm. However, for unconfined systems it is often 
useful to use a weak topology instead of the strong one. The trace-class 6 1 (S) k ) 
is the dual of the space JC($j k ) of compact operators, which is separable. Hence 
(3 1 (i^ fc ) can as well be endowed with the corresponding weak-* topology and its 
unit ball is sequentially compact for this topology. A natural question which 
arises in the case of a lack of compactness is then that of the weak-* limit of 

(k) 

the sets V N ■ We therefore introduce the set 

7>« := { 7 (2) € ©W), : 3 7$ € 7 « 7 «} . (L13 ) 

In Sectio n[2]w e will prove a weak version of the quantum de Finetti theorem (see 
Theorem 2.2) which implies that Vi, is the convex hull of all the weak limits 



of /c-particle density matrices of Hartree states and, therefore, any g < p^> 
can be written in the form 

7 (fc) = I (1.14) 
JBS) 

where [i is now a Borel probability measure on the unit ball BS) = {u £ : 
\\u\\ < 1} of the one-particle Hilbert space $j, instead of the unit sphere. 

Our Theorem |2.2| below is actually stronger and it states that if we have 
a sequence of iV-body states such that 7$ ^* 7^ for all k > 1, then the 
limiting density matrices 7^ can all be written as in (1.14) and they all share 



the same measure \i. By using this fact, we can easily prove the validity of 
Hartree's theory for systems in which the particles escaping to infinity carry a 
non- negative energy, that is, when the energy is a weakly lower semi-continuous 
function of the one- and two-particle density matrices. This is ensured if the 
particles may not form any bound states at infinity. We give interesting exam- 



ples of such systems in Section 4.2 below. In particular, we are able to provide 
a short proof for bosonic atoms, as considered first by Benguria and Lieb in [9]. 

The weak version of the quantum de Finetti theorem is related to recent 
results of Ammari and Nier [31 0] (see Section [2] for a more precise discussion) , 
but we follow a different approach, based on geometric methods for many-body 
systems [36] . Our proof of the weak quantum de Finetti theorem is in fact only 
the first step towards a more precise understanding of the lack of compactness 
for general systems and we will repeatedly use more refined arguments in the 
paper. 

Indeed, when the particles escaping to infinity have a nontrivial behavior, 
due for instance to attractive interaction potentials, looking at the weak limits 
of density matrices is not at all sufficient. The set V-w somehow only describe 
the particles which have not escaped, and the information on the other ones 
is completely lost. The accurate description of the lack of compactness will 
be done in this article, using the geometric methods of [36]. These couple the 
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(somehow algebraic) properties of many-particle systems with techniques from 
nonlinear analysis in the spirit of the concentration-compactness theory [39, 
147] . Our approach is now very general and it allows to cover many quantum 
systems, independently of the special form of their interaction. This is the main 
achievement of this article. For instance, we will recover the famous result of 
Lieb and Yau on boson stars [46J without using any particular property of the 
Newton potential. Note that our method is based on compactness arguments 
and it does not give any quantitative estimate on the discrepancy between the 
full many-body problem and its mean-field approximation, in contrast with 
operator-based methods that use specific properties of the interaction (see for 
example EE EU [57] ) . 

Typical result. For the convenience of the reader, we state now a typical result 
that can be obtained from our method, and which we will prove in Section [4] 
below. We consider the iV-body Hamiltonian 

N 1 

H v N =Y J {{m 2 -^ 3 ) s -m 2s + V{x 3 )) + w - l £ w(x k - x t ) (1.15) 

j=l l<k<£<N 

on the bosonic space Sj N = (g)f L 2 (R d ) with d > 1. Here m > and s £ (0, 1] 
are given constants. The case s = 1 corresponds to non-relativistic particles, 
whereas s = 1/2 describes a pseudo-relativistic system similar to the boson 
stars studied in [HI US]. To make a symmetric operator on $j N , as usual 
we assume that 

V : M d — > R and w : M d — > M are measurable and w(x) = w(—x). (1-16) 

The case of confined systems is the simplest, as explained above, and we refer 
to Section [3] where our results in this case are stated. Here we think of H% as 
describing an unconfined system so both V and w decay at infinity, which we 
formalize in the following 

Assumption on the decay at infinity of V and w. There exists some 
R > such that 

V1 RJ\B R =h + f2 and wl Kl! \ Bfl = h + U, (1.17) 

where fj G L Pj (M d ) with either max{l, d/(2s)} < pj < 00, or pj = 00 and 
fj — > at infinity, in the sense that |{s G R rf : |/j(^)| > &}| < 00 for all e > 0. 

We also need to specify which possible local singularities V and w may have, 
which is the subject of the following 

Assumption on the local singularities of V and w. There exists non- 
negative constants C , a± and (3±, with a_ + /3_ < 1, such that 

V±{x) < a±(-A) s / 2 + C and w ± (x) < /3±(-A) s / 2 + C, (1.18) 

where /+ = max{/, 0} and /_ = max{— /, 0} are respectively the positive and 
negative parts of f . 



Since, by (1.17), V and w are subcritical outside of the ball Br, the bounds 
in (1.18) are only interesting for the local parts V1b r and w1b r - Note that 
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under our assumption (1.18), the local singularities are allowed to be compa- 



rable to the kinetic energy and, in particular, V is not necessarily a compact 
perturbation of the kinetic operator. The upper bounds on V + and w + are not 
really necessary but they simplify the presentation. 

It is instructive to think of the caricature where both V and w are smooth 
functions of compact supports, which obviously satisfy our assumptions. The 
validity of Hartree's theory in this simple case is already a non trivial problem 
and does not seem to have been proven before. The conditions stated above 
are much more general however. In particular they are satisfied by Newton or 
Coulomb potentials when m > 0. 

Under the previous assumptions on V and w, the Hamiltonian is bounded 
from below and we denote by E v (JV) its ground state energy. The correspond- 
ing Hartree functional reads 



#(«) 



m 



A) 



a/2 



m s + Vj u 
1 

+ - 



w(x — y)\u{x)\ 2 \u(y)\ 2 dx dy 



and we denote by eg (A) its infimum on the sphere of radius y/X: 



inf £j[(u). 

ue.H"(R d ) 

\\u\\ 2 T o=\ 



Our main result contains two parts. The first item (i) deals with the validity 
of Hartree's theory at the level of the energy, independently of the strength 
of the external potential V (which may as well be V = 0). The second and 
third items of the statement give precisions about the density matrices of any 
sequence of approximate ground states, in particular when V is sufficiently 
negative to bind some (or all) of the particles. 

Theorem 1.1 (Validity of Hartree's theory). 



Assume thatV andw satisfy the previous assumptions (1.16), (1.17) and (1.18). 
(i) We always have 



lim 

7V-s>oo 



E V {N) 
N 



(1.19) 



(h) Denote by a sequence of approximate (normalized) ground states inS) N , 

that is, such that (^n^H^^n) = E v (N) + o(N), and byy^ the corresponding 
density matrices. Then there exists a subsequence (Nj)j>\ and a Borel proba- 
bility measure \i on the unit ball BSj = {u £ $) : \\u\\ < 1}, supported on the 
set 



M v = 
such that 



e£(l)-ek(l 



)}, (1-20) 



(h) 



\vP*)(u 9k \dn{u) 



M v 



(1.21) 



weakly-* in & 1 (Sj k ), for all k > 1. 
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(Hi) Assume now that the binding inequality 

e£(l)<e£(A) + 4(l 



A) 



(1.22) 



is satisfied for all < A < 1. Then the previous measure [i is supported on 

„(*) 



S$) and the limit (1.21) for ^y N is strong in the trace-class. In particular, if 

e^(l) admits a unique minimizer uh, u P to a phase, then there is complete 
Bose-Einstein condensation on it: 



7 <j=) _> \u% k ){u® k \ strongly m S 1 ^* 



(1.23) 



for any fixed k > 1 . 



Remark 1.2 (Generalizations). Our approach also applies for bosons in a mag- 
netic field corresponding to replacing the fractional Laplacian (m 2 — A) s by its 
magnetic version (m 2 + |V + iA(x)\ 2 ) s , see Remark 5.1 We only need that 



|^4| 2s satisfies similar assumptions as V. We are also able to deal with particles 



hoping on a lattice, see Remark 5.2 



It is a classical fact used in variational methods that the large inequality 

eg(l)<e£(A) + e&(l-A) 

is verified for all < A < 1. The set A4 V defined in (1.20) contains all the 
minimizers of the variational problems eK(A) for all A G [0, 1] which satisfy the 
equality 

e£(l) = e£(A) + e&(l-A). 

The interpretation of this condition is that a mass 1 — A can be sent to infinity 
without changing the lowest energy of the system. 



The role of the strict binding inequality (1.22) is precisely to prevent the 



particles to escape at infinity, and they are very often encountered in nonlinear 



models. Using Lions' terminology, the assumption (1.22) is here to avoid di 



chotomy, that is, to ensure that it is not favorable to split a minimizing sequence 
in pieces. 

In the translation-invariant case (V = 0), the many-body Hamiltonian does 
not have any ground state and there are sequences of approximate ground 

(k) 

states for which the density matrices j N all weakly tend to 0, even after a 
space translation (this is called vanishing in Lions' terminology). This is of 
course not in contradiction with (1.21) since in this case Ai 
u 



A4° contains 



0. On the other hand, there are other sequences (made of Hartree states 
for instance) which converge to a ground state of the Hartree functional, when 
it exists. 



It is interesting to note that the strict binding inequality (1.22) is only as 



sumed for the effective Hartree theory, and that it implies the expected behavior 



for the many-particle states. However, it is very important not to confuse ( 1.22 ) 



with the corresponding binding condition for the many-particle Hamiltonian 



N' 



By the HVZ theorem, the infimum of the essential spectrum of HY, is 



inf <j ess (iT^) 



inf 

fc=l,...JV 



(E v (N-k)+E°(k)), 



(1.24) 
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see, for example, [36, Thm. 12]. In particular, we have inf <r ess (iJ]Q < E v (N — 
1) and, applying the theorem to E v (N — 1), we find that 

as well. Therefore, the lowest eigenvalue E v (N) (when it exists) always behaves 
the same to the first order as the bottom of the essential spectrum. One has to 
go to the next order in the large- ./V expansion in order to distinguish the first 
eigenvalue from the bottom of the essential spectrum [37] . The Hartree binding 



condition (1.22) only counts whether it is interesting to send a number of order 
N of particles to infinity, whereas the HVZ criterion deals with any number of 
such particles, especially a number of order one. 

Note that the link between the Hartree minimizer (when it exists) and the 
many-particle ground states is only expressed here in terms of the density ma- 
trices ^y N . It is wrong in general that an approximate ground state ^jv £ Sj 
is close to a state u® N for the norm of Sj . One has to go to the next order in 
N to understand precisely the link between the two wave functions, see |37j . 



If there is complete Bose-Einstein condensation (BEC) as in (1.23), and if 
the unique Hartree ground state is non-degenerate, then it is shown by Lewin, 
Nam, Serfaty and Solovej in j37| that the energy can be expanded as 

E v (N) = Nel{l) + el + o{l) 

where eg is the ground state energy of an effective operator in Fock space called 
the Bogoliubov Hamiltonian. Since the validity of Hartree's theory and complete 
BEC were assumptions in [37], the present work supplements the article 1371. 

Before going to the more technical parts of the paper we summarize infor- 
mally its main message: 

(1) The validity of Hartree's theory in the mean-field limit may be viewed 
as a consequence of the structure of the set of bosonic iV-particle density 
matrices for large N. In this limit the representability problem can be 
given a satisfactory answer via the quantum de Finetti theorem. For 
confined systems, no other ingredient is needed. 

(2) For unconfined systems where particles are allowed to escape to infinity, 
a deeper analysis is required and it can be realized by combining a weak 
version of the quantum de Finetti theorem together with the geometric 
methods for many-body systems described in [36] . 

(3) Our method is general and can be applied in many different situations 
(non-relativistic or relativistic particles living in a domain, in the whole 
space or on a lattice, with or without external fields, etc). 

Organization of the paper. In the next section we discuss two versions of 
the quantum de Finetti theorem. After having recalled the usual statement of 
St0rmer and Hudson-Moody, we prove in Theorem 2.2| a weak version which 
will be very useful throughout the paper. In Section [3 we quickly explain how to 
deal with confined systems (in which there is no lack of compactness at infinity) , 
using the usual quantum de Finetti theorem. This is mainly introduced here 
for pedagogical purposes. Then, in Section[4]we study the Hamiltonian (1.15). 
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We start with the case of repulsive systems (for instance w > 0), for which the 
energy is a weakly lower semi-continuous function of the one- and two-particle 
density matrices. In this case the proof is an immediate consequence of the weak 
quantum de Finetti theorem. This covers bosonic atoms for example. Then, in 



Section 4.3, we investigate purely translation-invariant systems, by using some 
ideas of Lieb and Yau |46j . coupled to the geometric methods of [36J. Finally, 
we prove Theorem |1.1| in Section |4.4| and discuss some of its generalizations in 
Section [5j An alternative proof of the weak de Finetti theorem is presented in 
Appendix [Aj 

Acknowledgment. The authors acknowledge financial support from the Eu- 
ropean Research Council under the European Community's Seventh Framework 
Programme (FP7/2007-2013 Grant Agreement MNIQS 258023). 

2. A WEAK QUANTUM DE FlNETTI THEOREM AND GEOMETRIC 

LOCALIZATION 

The quantum de Finetti theorem is about the structure of bosonic states on 
the infinite tensor product of a given algebra. A simple formulation can be 
given in terms of density matrices only, following [31] . 



Theorem 2.1 (Quantum de Finetti). 

Let Sj be any separable Hilbert space and denote by Sj k := S) the correspond- 
ing bosonic k-particle space. Consider a hierarchy {7^}2Lo °f non - ne 9 a ti ve 
self-adjoint operators, where each 7^ acts on f) k . We assume that the hierar- 
chy is consistent in the sense that 

Tr fe+1 ^ +n7 ( fc+n )= 7 ( fc ) (2.1) 

for all k, n > 0. We also assume that = 1, which then implies Tr^k 7^ = 1 
for all k>0. 

Then there exists a Borel probability measure [i on the sphere SSj of $), in- 
variant under the group action of S 1 , such that 



7 (fc) = / \ u ®k)( u ®k\ d ^ 
JSS) 



(2.2) 



for all k>0. 

If moreover we have Tr(r7^ 1 )) < 00 for some self-adjoint operator T > on 
S), then n is supported in the quadratic form domain Q(T) of T . That is, we 
have fi(u) = for all u £ SSj \ Q(T). 

The result is the quantum equivalence of the famous Hewitt-Savage theorem 
for classical systems [HI [151 022 EH EH1 HE] ■ The latter deals with a hierarchy of 
symmetric probability measures fj,^ on £l N such that ^ k \A) = ^ k+n \Ax Q n ) 
for any k,n > and any measurable set A C 0, k . The quantum de Finetti 



Theorem 2.1 was proved in [65] [31] (see [261 02] for related content). The usual 
statement does not include the part concerning the operator T, but this part is 
easily shown by using the Hilbert space structure associated with the quadratic 
form of T instead of the original one. 
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If we are given a sequence of states f jv S Sj N with N — > oo, we obtain a 
sequence of density matrices, (7^ )o<k<N, which form a consistent hierarchy 
in the sense of (2.1 ), but only for k + n < N: 

\^n)(^n\, 




Tr fc+ i- 




(AO 



for < k < N, 
for k > N+ 1. 



(2.3) 



As each 7^/ is non- negative and has a trace normalized to 1, it is bounded in 
the trace-class &*(Sj ) for every fixed k > 1. Therefore it has a subsequence 
which converges weakly-* in the trace class, 7^ ^* 7^). By the diagonal 

procedure, we can make all the density matrices weakly-* converge along the 
same subsequence Nj — > 00. So we obtain in the limit an infinite hierarchy of 
density matrices (7^)fc>o, where each 7W acts on S) k . 

In general the sequence (^ k ')k>o is not consistent because the trace is not 
continuous for the weak-* topology in infinite dimension. Indeed, when passing 
to the weak limit we find by Fatou's lemma for trace-class operators 



7 ( fe ) = w-lim 7^ = w-lim Ti k+1 ^ k+n 7^ n) 
> Tr k+1 ^ k+n w-lim 7^~j~"' ) : 



Tr fc 



+l-S-fc+n7 



(k+n) 



where w-lim denotes the weak-* limit in Equality will not hold in 

general in the above equation and so we only have a priori 



y (fc) > Th 



7 VV > '^k+i^k+n ' "' , Vn, k > 0. (2.4) 

If Tr(7^ 1 ' 1 ) = 1 then it can be proved that Tr(7^^) = 1 for all k > 1 (see 
Corollary 2.4 below). The convergence must then be strong by the reciprocal 
of Fatou's lemma (see for example [60], Add. H]), and the sequence (7^ fc - ) )fc>o 
is consistent. In this case the quantum de Finetti Theorem 2.1 applies to the 
limiting hierarchy (7^' )&><)• 

In the systems in which particles are allowed to escape to infinity, a non- 
consistent hierarchy can be obtained in the limit. Consider, for instance, a 
given orthonormal basis (uj) of the one-particle space Sj, and the Hartree state 

= ( cos(#) ui + sin(0) uat)® 
for some 9 G [0, 2ir]. In this case we get 



(fc) 



cos(0) u\ + sin(6>) UN^ k \([ cos(6) u\ + sin(6>) 



cos(#) iii) 



cos 



(9) U1 y 



We see that in this example the sequence (7^) can be represented by a formula 
similar to ( |2.2[ ), but with a measure [i that is a 5 at the point cos(#) u\ in the 
unit ball BS) = {u G fj : \\u\\ < 1} of $). That we end up with a measure 
living on the unit ball BS) instead of the unit sphere SSj is not a big surprise, 
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of course, as we are considering weak limits. This is actually the general case, 
as stated in the following result, which will be an important tool in this article. 

Theorem 2.2 (Weak quantum de Finetti). 

Let Sj be any separable Hilbert space and denote by $j k := ® s $j the correspond- 
ing bosonic k-particle space. Let be any sequence of mixed states on Sj N 
(that is, Ttv > and Tr^jv Tj^ = 1) such that 

v( fc ) 



(k) 



r 



weakly-* in the trace class © 1 (Sj k ) for all k > 1. Then there exists a Borel 
probability measure [i on the unit ball BS) of ft, invariant under the group 
action of S 1 , such that 



dfi(u) 



u 



><«« 



(2.5) 



for all k > 0. 

Moreover, if we have Tr (Tj^) < oo for some self-adjoint operator T > 
on Sj, then /j, is supported in the quadratic form domain Q(T) of T, namely 
H(u) = for all u£ BS)\ Q(T). 

Ammari and Nier have recently proved in [31 0] results that imply Theo- 
rem 2.2 In analogy with semi-classical analysis, they called fi a Wigner mea- 



sure. They deal with an arbitrary sequence of states in Fock space and, there- 
fore, obtain in the limit a measure fi which can live over the whole one-particle 
Hilbert space fj, instead of the unit ball as in our situation. In [3l Thm. 6.2] 
they first construct the Wigner measure /i by testing against anti-Wick and 
Weyl observables, before looking at Wick observables which are related to den- 
sity matrices in [31 Cor. 6.14]. The case of [i having its support in a ball 
is studied in [1]. The connection between Wigner and de Finetti measures is 
discussed in [2, Section 6.3]. 

In the present paper we provide two different proofs of Theorem|2.2[ which are 



both based on Theorem 2.1 The first proof is based on the finite- dimensional 
de Finetti Theorem and on the geometric techniques introduced in |36j and it 
has the merit of clarifying how the measure /i arises in case the density matrices 

(k) 

7^ do not converge strongly. This will be particularly important to understand 
unconfined quantum systems in the rest of the paper and hence we explain this 



first approach in details in this section. It is also possible to prove Theorem 2.2 



by following arguments similar to those of Hudson and Moody [31], and we 
quickly explain this in Appendix [A] for completeness. 

Remark 2.3. It could seem an interesting question to find the structure of the 



set of all the density matrices satisfying the inequality ( |2.4[ ). In Theorem 2.2 
we only characterize those arising from a sequence of states (r^r) with - 
The set of density matrices satisfying (2.4) is way too large, however. 



oo. 
For 



instance it contains the density matrices of n-particle states with n fixed, for 



which a representation of the form (2.5) cannot hold. Consider 

\u)(u\, 7 (n) = forn>2, 



7 



(0) 



1, 



7 



(i) 
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which of course satisfies ( 2.4 ) . These density matrices are those of a one-particle 



state u £ Sj and they cannot be written in the form (2.5). The measure [i would 
need to be ^ = 5 U but then ^ for n > 2. 



By using the weak de Finetti Theorem 2.2 we easily recover the well-known 



(k) 

fact that the density matrices 7p^ all converge strongly if and only if the one- 
particle density matrix 7^ converges strongly. 

Corollary 2.4 (Strong convergence). 



Let {iV} be as in Theorem 2.2. Then the following are equivalent: 

(1) Tr fi 7 « = 1; 

(2) 7^ — > 7W strongly in the trace-class (S 1 ^); 

(3) Tr fifc 7W = 1 for all k > 1; 

(4) 7^ — > strongly in the trace-class & 1 (Sj k ) for all k > 1; 

(5) the measure of Theorem \2.2\ has its support on the sphere SSj. 

Proof. The equivalence between (1) and (2) (and (3) and (4)) is the reciprocal 
of Fatou's lemma mentioned above (see, for instance, [15] , [Ml Cor. 1] and [6"0"1 



Add. H]). The rest follows directly from (2.5). □ 



Geometric localization. In this section we provide a proof of Theorem 2.2 
using the concept of geometric localization^ in Fock space, which was developed 
in [17\ Q] and thoroughly used for nonlinear many-body systems in |36| . This 
method is the first step to understand the lack of compactness of the sequence 
(&n) °f approximate minimizers appearing in Theorem 

For any mixed state Tat in the iV-particle space Sj N and any localizing op- 
erator < A < 1 on Sj, the corresponding localized state is by definition the 
unique state in Fock space 

F{Sj) = c e ^ e s) 2 © • • • 

which has -particle density matrices, for all 1 < n < N. The 

localized state lives in the truncated Fock space C © Sj © • • • © $3 © • • • and 
it is often a mixed state, even when the initial state r^r is pure. This localized 
state can be explicitly computed (see [36J, Example 10]) and it is equal to 

Gn = G Nfl © G N1 © • • • © G NjN © © • • • (2.6) 

where 

a ( N\ / ^,7 / ®N-k „, / ®7V-fc 

G^ k =() Tr k+1 ^ N (A 9k © T N A® k ® 



An important property of the localization of iV-particle states is that 

T N,N-k 



(2.7) 



l^fi* G N b — Tr^iv-fc GJ~ N \ (2. 



^The term "geometric" was first used in [591 158] to denote the use of partitions of unity in 
the configuration space for many-body systems. 



16 



M. LEWIN, P.T. NAM, AND N. ROUGERIE 



for all k = 0,...,N. If A = Id for some D C M. d , then (2.8) simply means 
that the probability of having k particles inside D is equal to the probability of 
having N — k particles outside D. 
A simple calculation also shows that 



N 



n 



-1 N 

£ 

k=n 



n 



Tr r , 



(2.9) 



see [36] for details. We emphasize that in [36; a different convention is used for 
the n-particle density matrix, and this is responsible for the additional factor 

Remark 2.5. When A = P is an orthogonal projection, geometric localization 
is the same as using the isometry of Fock spaces J~(S)) = J- (PS)) J-(P^S)) 
and restricting the state to the smaller Fock space J- (PS)), by taking a partial 
trace in the second variable. 

The link between the de Finetti measure /i and geometric localization is 
emphasized in the following result, which will be very useful in the proof of 
Theorem 11.11 below. 



Theorem 2.6 (De Finet ti m easure and geometric localization). 

Let (r^v) be as in Theorem 2.2 and assume, furthermore, that TrT7p^ < C for 
some non-negative self-adjoint operator T. Let < A < 1 be an operator on S) 
such that A(l + T) -1 / 2 is compact and let G^^be the associated localized state 



in the Fock space F(S)) defined in (2.6) and (2.7). Then 

f (£) Tr fife G A ^ k = I d^(u) f(\\Auf 



lim 



k=0 



(2.10) 



for all continuous functions f on [0, 1]. 



Here geometric localization is used to detect the particles which do not escape 
to infinity and the role of the operator A is to turn weak convergence into strong 
convergence |36j . The interpretation of the convergence (2.10) is that the mass 
of the de Finetti measure [i on the sphere {||u|| 2 = A} is the probability that a 
fraction A of the particles does not escape to infinity. 

Of course, we can turn the matter around and obtain an information on the 



number of pa rticles that have escaped. Using the fundamental relation (|2.8|), 
Theorem 2.6 implied 

N 



lim 

N-+oo 



k=o v 7 



Tr^fe G 



N.k 



lim 



N 

X 

k=0 



/ I 



A: 
N 



Tr fifc G A 



N.k 



dfi(u) f(l - \\Au\\ 2 ). 



(2.11) 



(1 



°The convergence (2.101 does not hold in general if A is replaced by (1 
- A 2 ) 1/2 (l + T) -1 ' 2 
< I- 



is not always compact. Indeed, 1 — ||^t|| 2 > || Vl — A^i 



A 2 ) 1/2 , because 



when 
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In practice, we will use (2.11) with A = xr{ x )i a localization function in a ball 



of radius R, which is compact relatively to the fractional Laplacian. 



The rest of this section will be devoted to the proof of Theorem 2.2 and 



Theorem 2.6 using geometric localization. 



Proof of Theorem 2.2. Our approach can be summarized as follows. First we 
localize the quantum state in a finite-dimensional space using a finite dimen- 
sional projection P, in order to convert the weak convergence into the strong 
convergence. By doing so we get the localized state in Fock space and we 
apply the quantum de Finetti to each of the projections k in the fc-particle 
spaces. The limiting density matrices have a representation in "spherical coor- 
dinates" on BS"), 



7it )= / / dv loc (\,u) \ n \u® n )(u® n \, (2.12) 

JO JSS) 

where A corresponds to all the possible values of k/N. In the end of the proof 
we remove the localization and get the result. 

To carry out this program, we fix an orthogonal projection of finite rank P on 
fj and let be the corresponding localized state as defined in (2.6) and (2.7). 



We write the sum in (2.9 ) as an integral over an additional parameter < A < 1, 
which is A = k/N. Let MpL be the following Radon measure on [0, 1], with 



values in the set of self-adjoint operators in @ 1 ( 
matrices of size dim <g>™(Pi})): 

N 

r(n) 



*P$j) (that is, hermitian 



dM^(A) := ^ 5 k/N (X) Tr n+1 ^ k G 



N,k- 



This measure satisfies 
fi 



/ 

Jo 



N 



N 



TV dil4 n i(A) 



G 



< 



G 



p 

N,k 



Then, using (2.9), we have for all n > 



Tr 



T jv 



\ n dM^ N (\) 



k=0 



N 



k=n 



We can write 

k 
N 



N 



n 



and by using Bernoulli's inequality 



N 



n 



-i 



n 



j_ 

N 



2_ 

N 



dM^ N (X) = 1. 



TrG p 



N,k- 



N 
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we obtain 
Therefore, 



-i 



< 



N 



N 



in 



1) 



< 



in 



If 



N -n+l ~ N - n + 1 



(2.13) 



Tr 



p®n An) p®n 



A n dM (n) 



(A) 



in 



1) 



k=n 



Since P is a projection of finite rank and by assumption 7^ — weakly-* in 
the trace-class, the sequence p® n ^^p® n converges strongly in the trace-class 
^ p®n^(n) p®n _ Therefore we have proved that, for any fixed n, 

■1 



N — n + l 
(2.14) 



lim Tr 







A»dil4 n i(A) 



0. 



By construction, M p n ^ N is a bounded sequence of (Radon) measures on [0, 1] with 
values in the cone of non-negative operators in the finite-dimensional space 
(&™PS) C Sj n having a trace < 1. This sequence has a subsequence which 
converges weakly (in the sense of bounded measures on a compact set) to some 



(n) 

measure M p and is therefore such that 



X n dM { p n) (X) 



for all n > 1. For n = 0, M { p ] is a Borel probability measure on [0,1]. For 

(n) 

n > 1, the value of the measure M p at A = is not important, as it does not 
contribute to the density matrices. In order to simplify our reasoning, we will 
simply take M p n) ({0}) = M p 0) ({0}) Id 8 « Pj5 Tr (ld^p^)~\ 



(n) 

Now, we claim that the sequence M p is consistent in the sense that 



T^n+l-tn+k M 



(n+fc) 
P 



M 



(n) 



(2.15) 



for n, k > 0. Indeed, we have 



Tr 



n+l— >n+k Mp jy 



(n+k) 



N 



(A) 



8j/nW Tr„ + i^j Gat,. 

j'=n+fc 



4(A) 



n+fc— 1 

E 

j=n 



and therefore 



A Tr c 



Tr n+1 ^ n+fe M}T+ fc) (A) - M^(A) 



(n) 



. , -. n+fc— 1 

n + fc — 1 



^■ /7V (A) TrG p 



N,j- 



j=n 



The sum of the right side is a uniformly bounded measure on [0, 1], hence we 
obtain fl2.15h. 
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By the (finite-dimensional) quantum de Finetti Theorem 2.1 there exists a 
Borel probability measure vp on [0, 1] x SSj n (PS)), invariant under the action 
of S\ such that 

dM { p n) (X)= [ dv P (\,u)\u® n )(u® n \. 
JSS) 

The original statement applies only for each fixed A G [0, 1]. In order to deal 
with the present case of a probability measure on [0,1], we first approximate 
Mp (A) by a step function and then we pass to the limit. The weak limit is 
still a Borel probability measure as we are working on a compact set of a finite 
dimensional space. We have now proved that 

p®n 7 (n)p®n = f f ^ P (A, «) A n |u®™) 
JO JsS) 

= 11 dv P {\,u)\(V\u)® n )((V\u)® n \. 
Jo Jssj 

Associated to the probability measure up on [0, 1] x Sfy n Pfj, there is a unique 
probability measure /ip on B$j n PSj which is formally given by the formula 
Hp(u) := vp(\\u\\ , u [|it||~ ). More precisely, fip is defined by its action on 
continuous functions on Bfj by 

f(u)dpL P (u) := / / f(yf\u)dvp(\,u) 



IBS) Jo JSS) 

and we get as we wanted 

for all n > 1. 

The argument can be applied for any chosen finite-rank orthogonal projection 
P. Consider now a sequence Pk which converges strongly to the identity on $j 
and apply the previous argument for each k. We find a sequence of Borel 
probability measures /i& on BP^Sj, invariant under the action of S , such that 

p®n 1 (n)p®n = f dfJL k (u) \u® n ) {u® n \. 
JbS) 

By construction the measure coincides with \m for I < k on cylindrical Borel 
sets having their base in Pi$). Since these measures all have their support in 
a bounded set in S), there exists by |61l Lemma 1] a unique Borel probability 
measure \i on Bf) which coincides with on cylindrical Borel subsets having 
their base in -Pfcfj. This precisely means that 

d» k (u)\u® n )( u ® n \= [ d^u)\{p k ur n ){{p kU r n \ 

BS) JBS) 

and therefore 

pm^pm = p®n (^j d ^ | u ®n^ u ®n|^ p»n_ 

Taking now fc — )• oo finishes the proof. □ 
We now turn to the 
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Proof of Theorem 2.6. It suffices to prove the result for /(A) = X p with p = 
0, 1, .... For p = 0, this is nothing but the fact that X^fcLo r ^ T ^ k ^Nk = 1 since 
G^t is a state. For p = 1, we note that 



r 7V 

AT 



A; 



£^Tr^G£ jfc = Tr(A 7 gA) 



fe=0 



as follows from (2.9). The right side can be written 

Tr(A 7 ^A) = Tr(A{T + 1)~V 2 (T + l) 1 /^^ + l) 1 / 2 ^ + l)- 1 / 2 ^) 

and this converges to Tr(Aj^ A) since A(T + l)^ 1 / 2 is compact and (T + 
1)V 2 7 W (T + 1) 1/2 ^* (T + l) 1/2 7 (1) ( r + 1) 1/2 weakly-* in 6 1 . So we obtain 

^°°fc=0 



d/i(u) Pn|| 2 . 



The proof is similar for p > 2, using that 



and the estimate (2.13). □ 



3. Validity of Hartree's theory for trapped bosons 



As a direct application of the usual de Finetti Theorem 2.1, we state here a 
rather general result for trapped systems, which is in the spirit of the earlier 
works [231 [501 [52]. 

Let $) be a separable Hilbert space and consider an iV-body Hamiltonian of 
the form 

N 1 

j=l l<k<£<N 

acting on the bosonic space Sj N = t^)^ fj . We make the assumption that 

T is bounded from below and it has a compact resolvent (3-2) 

which is the mathematical formulation of the system being trapped. We also 
assume that w is a self-adjoint operator on i^ 2 which is relatively small compared 
to the one-particle term, in the sense that 

- /3_ (T ® 1 + 1 <g> T) - C < w < p + (T ® 1 + 1 <g> T) + C (3.3) 



where C > 0, /3+ > and 1 > /3_ > 0. The upper bound in (3.3) is not essential 
but it simplifies the analysis. Also, the precise form which we have chosen for 
Hn is not really important. We could as well introduce 3-body terms, or use 
abstract conditions like the ones in 1:521 [701. 
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Under the assumption (3.3), we deduce that 

N 



A' 



(1 - /3_) Tj -CN<H N <(l + p)Y^ Tj + CN. (3.4) 
j=i 3=1 

This proves that Hn is bounded from below and that its Friedrichs extension 
has the same form domain as the one-particle term ^2f = i Tj. 

3.1. Zero temperature case. At zero temperature, we are interested in the 
ground state energy of Hn, which is given by 



E(N) = mia^ N (H N ). 



(3.5) 



The assumption that T has a compact resolvent easily implies with (3.4) that 
the spectrum of Hn is purely discrete, hence E(N) is an eigenvalue of finite 
multiplicity. 

The corresponding Hartree problem is given by 



e H := 



ueQ(T) 
||«||=i 



ueQ(T) 

\H=i 



inf £n(u) := inf < (u, Tu) + - (u ® u, w u ® u) s 



(3.6) 



where Q(T) is the quadratic form domain of the bounded-below operator T. 
Again the compactness of the resolvent of T and our assumptions on w easily 
imply the existence of at least one minimizer for en- Also, the set of minimizers 
M C Sf) is bounded in Q(T). 

Theorem 3.1 (Validity of Hartree and BEC for trapped bosons). 

Under the previous assumptions (3.2) and (|3.3l) on T and w, we have 



Km ^ 



e H . 



If O^iv) is any sequence such that {^n, Hn^n) = E(N) + o(N), then there 
exists a subsequence and a probability measure fi on the set A4 of minimizers 
of en (modulo a phase), such that 



1 i m t£. 



®&\ /„,®fcl 



M 



strongly in the trace- class for any fixed k and, even, for the norm induced by 
the quadratic form ofT: 



lim 



A: 



1/2 ( (fc) 



M 



dniu) \u® k )(u® k \ J 



0. 



(3.7) 



where A^ := X^=iC^ + C). in particular, if en admits a unique minimizer uq, 
then there is complete Bose-Einstein condensation on uq: 



lim 7 £ = !<*>«' 



Mk\ /„,®fc| 



(3- 



Theorem 3.1 applies to several typical situations. 
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Example 3.2 (Bounded domain). Take f} = L 2 (f2) where is a bounded 
domain of M. d and T = —A + V with chosen boundary conditions and w^e = 
w(xk — xi) with V, w G L p (£l) for some p > max(l,d/2). Hartree's theory is 
then valid in the mean-field limit by Theorem |3.1[ 

In one dimension, the interaction can also be a delta potential, w^g = \5(xk— 
Xi), as in the Lieb-Liniger model [30] on a finite interval Q = (0,1). This 
potential, which acts as 



0M(xi -aj 2 )¥) := / |*(s,s)rds (3.9) 
Jo, 

on two-body wave functions \&, is relatively form-bounded with respect to the 



Laplacian similarly as in (3.3), by the Sobolev embedding. By Theorem 3.1 



the associated Hartree functional 

£ u (u) = (\u'{x)\ 2 + V{x)\u{x)\ 2 + ^\u{x)\^j dx 

gives the leading order of the energy per particle when N — > oo, for any fixed 
A G M. For A > this has been proved before in [57] . 

Example 3.3 (Confining potential). Take S) = L 2 (R d ), T = -A + V(x) 
with V G L^ oc (M. d ) satisfying hni|. c |_ 5>00 V(x) = +oo, and w^t = w(xk — xg) with 
w G Lf oc (M d ) satisfying 

-a(V(x) + V(y)) - C < w(x - y)t(\x — y\ > R) < P{V(x) + V{y)) + C 

for some constants p > max(l, d/2), C > 0, R > 0, j3 > 0, 1 > a > 0. We see 



that the above assumptions <\3.2h and (3.3) are satisfied, and therefore Hartree's 



theory is again valid by Theorem 3.1 



When w > 0, this was proved in [2 



where the energy E(N) is even expanded to the next order. 

Example 3.4 (Magnetic fields, rotation). Take f) = L 2 (R d ), T = —(V + 
iA(x)) 2 + V(x) with V and w as in the previous example, and A G L^ oc (M. d ). 



Hartree's theory is again valid in this case by Theorem 3.1 One may for 
example consider A{x) = (B/2)(x2, — x±,0) the vector potential corresponding 
to a uniform magnetic field curLA of strength B pointing in the x% direction. 
Using the analogy between the Coriolis and the Lorentz force, one may also 
think of A(x) = (Q/2)(x2, —xi, 0) with V = V — \\A\ 2 , which corresponds to a 
gas rotating at speed around the X3 axis, in the trapping potential V . The gas 
is then described in the rotating frame and the potential V (called the effective 
potential) takes into account the effect of the centrifugal force. 

Example 3.5 (Bosons in the lowest Landau level). It is possible to deal 
with 3D bosons rotating along a fixed axis, in a harmonic confining potential 



and with a delta interaction (defined similarly as in (3.9)), provided that the 
Hamiltonian is restricted to the lowest Landau level. See for instance [44] and 
the references therein for a precise definition of the model. In the Bargmann 
representation, the one-particle operator becomes T = zd z where z G C is the 
complex coordinate in the plane orthogonal to the axis of rotation. The delta 



potential is itself a bounded operator on the Bargmann space. Theorem 3.1 
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applies to both the repulsive and attractive cases. The former was treated 
in 



We can now provide the 



Proof of Theorem\3l\ Let (V N ) be such that (^ N ,H N ^ N ) = E(N) + o(N). 

Up to a subsequence we may assume that 7^ — 1 7^ weakly-* in the trace class 

as N — > 00 for all k > 1. By (3.4) we know that Tr(T7^) is bounded and we 

therefore also get that (T + Cq) 1 / 2 ^ (T + Q)) 1/2 (T + C ) 1/2 ~f W {T + C ) 1/2 
weakly-* in the trace class, where Co is chosen such that T + Cq > 1. From the 
compactness of (T + Cq) -1 it follows that 



Tr 7 « =Tr(T + C r L (T + Q 



-1 



> 1/2 7^ 



(T + C ) 1/2 ->Tr 7 «. 



By Corollary 



2.4 



we infer that 7^ — > 7^) strongly in S 1 ^^) for all > 1. 
By Theorem 2.1 there exists a probability measure fi on Q(T) Pi Sio such that 



7 



(A:) 



®fc\ /„,®fc| 



for all fc > 1. Now we go back to the iV-body energy and notice that 
(* N ,H N * N ) _^ ( (I)' 



TV 



(a - 2)C + (1 - a) Tr fl (T + C )7* 



H — Tr fi 2 
2 fl 



(w + a(T <g> 1 + 1 <g> T) + 2C ) 7* 



(2) 



(3.10) 



where the constant Co is now chosen such that T + Co > 1 and w + a(T 1 + 
1 (g> T) + 2Cq > 0. By Fatou's lemma, we obtain 



lim 

N -5- 00 



N 



1 



v( 2 ) 



SSj 



d(i(u) £k(u) > e H . 



Recalling the upper bound E(N) < Neu, we conclude that 

dfi{u) £r{u) = e H . 



lim 

jV-s-oo 



N 



It is now clear that \i must have its support on the set M containing all the 
minimizers of en. Going back to the proof we see that equality must hold 
everywhere, which implies in particular that Tr(T7^ ) — > Tr(T7^ 1 )). By the 
reciprocal of Fatou's lemma, this gives the strong convergence 

(r + Co) 1 / 2 ( 7 S-7 (1) )(r + c ) 1 / 2 ^o 

in the trace class. To get the same for the higher density matrices, we use the 
same argument and the fact that Tr^ (T'y^J = k^ 1 Tr^k ^( Yli=i ^7*jv)' 

We mention that everything holds here for a subsequence. However, if Ai = 
{uq} then our proof shows that all the subsequences must have the same limit 
and the corresponding convergence follows for the whole sequence 7^, . □ 
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3.2. Positive temperature case. Our result can be extended to the case of 
positive temperature and we quickly explain this now. We assume that 

Tre^ T <oo (3.11) 

for all P > (which implies that the resolvent (T + C)~ l must be trace-class, 
hence compact). The condition can be relaxed by assuming only Tre~@ T < oo 



for one (3 > 0, but we use the stronger assumption (3.11) for simplicity. 



The iV-particle quantum free energy at temperature j3 1 > is now given 
by the formula 

1 



E(P,N) 



P 



log Tr 



(3.12) 



From (3.4) we have 



Tr. 



which proves that 



(l-a)T 



N 



1 



E(P, N)>-N[C+^log Tr^ e^ 1 "^ 



(3.13) 



We see that, under the assumptions (3.11) on T and (3.3) on w, the iV-particle 



canonical free energy E(f3, N) is well defined for all (5 > 0. Let us denote by 

P -I3H N 



0,N 



Tr fi jv e-P H * 



(k) 

the associated Gibbs state and by 7^ N = Tr^+i- 
matrices. Our result is then the following 



(3.14) 

>N^p,N its fc-particle density 



Theorem 3.6 (Trapped bosons at positive temperature) 



Under the previous assumptions (3.11) and (3.3) on T and w, we have 



lim 

N^foo 



E{P,N) 



N 



en 



for any fixed (3 > 0. The other conclusions of Theorem 3.1 remain valid for the 



k-particle density matrices 7g jy- 

Remark that, at a fixed temperature > 0, almost of the particles still 
condense on the Hartree ground state. In this case, the effect of the temper- 
ature is only seen on the next order of the free energy, given by Bogoliubov's 
theory (see, for example, [66], Prop. 3.1] for non-interacting systems and [37] 
for interacting systems). One can see the effect of the temperature in Hartree 
theory by taking a temperature which diverges proportionally to N, that is 
Pn ~ 1/N, see 



Proof. The proof is very similar to the one of Theorem 3.1 and we only explain 
the differences. By Gibbs' variational principle, we have 



E(P,N) 



inf (Tr^iv H N T + /T 1 Tr^iv T log r) 



TiT=l 
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Let r^ jv be as in (3.14). By arguing similarly as for the proof of (3.13), it is 
easy to verify that 

for a constant which is independent of N. Using this bound we can write 
TV H N T/3 tN + /3" 1 Tr^iv Tp^ logTp >N 

> E(N) + e TV \ J2 r /3,^v + /3" 1 TV T P,N log Tp, N - eCN. 

where we have used that Hn > E(N), the ground state energy at zero temper- 
ature. In order to control the error terms, we are now going to use the exact 
behavior of the free energy in the non-interacting case. 

Lemma 3.7 (A uniform lower bound in the non-interacting case). 

Let K be an operator on S) such that Tr^ e~@ K < oo and let Ai be the first 
eigenvalue of K which we assume to be non- degenerate. Then we have 

< -log (TV e-££i*) - NM{K) + log (xr, ± e _ (i J Al) _ J ^ 

(3.15) 



where S)± = l( Xl>00 )(K)S) C ft. 

The proof of the lemma is a well-known simple calculation which can be 
found, for instance, in |66|, Prop. 3.1]. There is a similar statement (with a 
different lower bound) when the first eigenvalue is degenerate, which we will 
not use here. 

In our case we do not know the multiplicity of the first eigenvalue of T, but 
we can argue as follows. First we add a positive finite-rank operator B to T in 
order to remove the degeneracy of \±(T), without changing the first eigenvalue. 
We obtain an operator T' = T + B to which we can apply Lemma 3.7 The 
error made by replacing T by T" is e \\B\\ N, leading to the following estimate: 

n r^jv log r^jv 



Now we divide by N and let N — > oo. By Theorem 3.1 we get in the limit 
liminf E i< 3 ^l > eH + £ (X 1 (T) -C- \\B\\). 



Taking e — > gives the result. The rest is proved exactly like in Theorem 3.1 

□ 
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4. Validity of Hartree's theory for unconfined bosons 

In this section we deal with quantum systems in which particles are allowed 
to escape to infinity. We will prove Theorem |1.1| which was stated in the 
introduction. 

4.1. Preliminaries. We consider the Hamiltonian 

N 

^:=E((™ 2 - A ;) S -™ 2S + ^)) + ^3T £ 

3=1 l<k<e<N 

on the iV-particle bosonic space fi N = (g)f L 2 (R d ), where m > and s G (0, 1] 
are given constants. It will be convenient to introduce the notation 

K = (m 2 - A) s - m 2s and T = (m 2 - A) s - m 2s + V. 



For the rest of the paper, we always work with Assumptions (1.16)— (1.17)- 



(1.18) on V and w. Under these assumptions, we have 

N N 

(1 - a_ - P-) Kj -CN <Hl <{l + a + + /3+) ^ Kj + CN (4.2) 
3=1 i=i 



for some constant C which could be different from the one of (1.18). In par- 
ticular, the quadratic form associated with HYt has the same domain as that 
of the free kinetic energy 

Yjj=i Kj, We denote by E V (N) the bottom of the 
spectrum of H^- The following says that E v (N)N^ 1 always has a limit. 

Lemma 4.1 (Monotonicity of the energy per particle). 

The sequence N i— > E v (N)N^ 1 is increasing and < 0. 

Proof. We have already explained in the introduction that the energy per par- 

(k) 

ticle is increasing because of the monotonicity of the sets V N . It remains to 
show that E V (N) < for all N > 1. By taking a Hartree state u® N with 
u e H s (R d ), we find 

E V (N) . „ . 1 

We fix a smooth normalized function of compact support v G C^°(M rf ) and 
take n = k^ 1 / 2 ^2j = i v{- — Xn). In the limit where \x 3 n \ — > oo as n — > oo and 

|a^. — x* n | — >• oo for j ^ j' , the energy of u becomes the sum of the energies of 
the pieces, leading to 

^jp- < ( v ' Kv ) + ^ ff w(x-y)\v(x)\ 2 \v(y)\ 2 dxdy. 
We have used here that 



< (u,Tu) + - // w(x — y)\u(x)\ 2 \u(y)\ 2 dx dy. 



lim // w(x — y)\v(x)\ \v(y + z)\ dxdy = 
1*1 ->•<». 



for any fixed function v G C^°(R d ), which follows from our assumption (1.17) 
on w at infinity. Taking now k -> oo leads to E V (N)/N < inf a((m 2 - A) 8 / 2 - 
m s ) = 0, which concludes the proof. □ 
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Now we turn to the elementary properties of Hartree functional 
S&(u) = (u, ((m 2 - A) 3 / 2 -m s + V^u^ + ^JJw(x- y)\u(x)\ 2 \u{y)\ 2 dxdy. 

By ( |4.2[ ), it is bounded from below on bounded subsets of L 2 (M. d ). We denote 
its infimum on the sphere of radius vA by 



e£(A) := inf f£(u). 

||«|| 2 =A 

Some basic properties of Hartree's theory are given in the following lemma. 

Lemma 4.2 (Hartree theory). 

Under Assumptions (1.16), (1.17) and ( 1.18[ ), we have 

e£(l) < e£(A) + e°(l - A) < e° H (A) < (4.3) 

/or all A G [0, 1]. Moreover, if 

e £(l)<e£(A) + e£(l-A) 

/or a// A G [0, 1), i/ien aZZ i/ie minimizing sequences for en(l) are relatively com- 
pact in H s (R d ) and converge, after extraction of a subsequence, to a minimizer 
inSSjf]Q(T). 



4.1 



Proof That eg (A) < e&(A) < was already shown in the proof of Lemma 
The rest follows from the concentration-compactness method and the localiza- 
tion procedure explained in the appendix of [35]. The proof will be omitted. □ 

In the next subsections we prove the validity of Hartree's theory for the many- 
particle operator HYr. We first deal with repulsive systems, then consider the 
purely translation-invariant case V = 0, before we turn to the general case. 

4.2. Weakly lower semi-continuous systems. We deal here with the case 
in which the particles escaping to infinity cannot form bound states. This is 
formalized by saying that E°(N) = for all N > 2 which, in the mean-field 
regime studied in this paper, turns out to be equivalent to E°(2) = 0. 

Lemma 4.3 (Absence of bound states at infinity). 

The following are equivalent: 

(i) the operator (m 2 — A) s — m 2s + w has no negative eigenvalue on L 2 (K d ); 
(ii)E°(2) = 0; 

{Hi) E°(N) = for all N > 1. 

Proof. Since E (2) is the bottom of the spectrum of H®, E°(2) = is clearly 
equivalent to = (m 2 — A x ) s + (m 2 — A y ) s — 2m 2s + w(x — y) > 0. By 
removing the center of mass similarly as in |38|, Appendix A], this is the same 
as [m — A) s - m 2s + w > 0. Finally, since iV i— > E°(N)N 1 is non-decreasing 



and < by Lemma 4.1 E°(2) = clearly implies E°(N) = for all N > 2. □ 



Now, the following says that, under one of the equivalent assumptions of 
Lemma |4.3[ the energy is a weakly lower semi-continuous function of the one- 
and two-particle density matrices. 
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Proposition 4.4 (Weak lower semi-continuity of the energy). 

We assume that V and w satisfy (1.16)-(1-17)-(1.18), and that w satisfies one 



of the equivalent assumptions of Lemma 4-3 Then the energy of H 



is a 



weakly lower semi- continuous function of the first two density matrices: If we 
have two sequences 7]^ , 7^ > with Tr^2 7]^ = 1 and 7*- 1 ** 
* 7W weakly-* in & 1 (Sj k ) for k = 1,2, then 



that 7^ 



Tr2 7jv > suc h 



liminf ( TreJT7 



+ 2 Tl >> 2 Kx 



|) >Tr fi [T 7 «] + ^Tr 



Sj2 [wy 



(2)1 



(4.4) 



By using Proposition 4.4 we can now easily prove the main theorem for such 
systems. 

Theorem 4.5 (Validity of Hartree's theory in the wise case). 

We assume that V and w satisfy (1.16)-(1.17)-(1.18), and that w satisfies one 



of the equivalent assumptions of Lemma 4-3 
(i) We always have 



lim 

TV-s-oo 



E V (N) 
N 



(4.5) 



(n) Denote by a sequence of (normalized) approximate ground states inSj N , 

that is, such that (*jv? -ffjv^iv) = E v (N) + o(N), and byyjp the corresponding 
density matrices. Then there exists a subsequence (Nj) and a Borel probability 
measure \x on the unit ball BSj, supported on the set A4 V = {u £ BS) : £^(u) = 
e H(IM| 2 ) = e H(-0}; such that 



7jv. 



\u® k ){u® k \dn(u) 



(4.6) 



weakly-* in & l (Sj k ), for all k > 1. 



(Hi) Assume now that the binding inequality e^(l) < e^(A) is satisfied for all 
< A < 1. Then the previous measure \i is supported on S$) and the limit (4.6) 

(k) 

for j N is strong in the trace-class. 



The proof of this theorem is an immediate consequence of Proposition 4.4 
and of the weak de Finetti Theorem 12.21 



Proof. Let {^tv} be a sequence of approximate ground states as in the state- 
ment (ii). Extracting subsequences if necessary, we may assume that 7^ 
7W for all k > 1. By the weak de Finetti Theorem 
probability measure fi on BSj n H s (R d ) such that 



2.2 



there exists a Borel 



7 



(k) 



BSj 



d/i(u) |u® fe )(u® fe | 



(4.7) 
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for all A: > 1. From Proposition 4.4 we have 



lim ^> 



lim ( Tr fl [T>y$] + Tr^ [wj 



(2)- 

N - 



>T^[T^] + Tv^[w^]= [ dn{u)£V(u)> [ dn{u)el{\\u 



By (4.3) we have epj(||n|| ) > ejj(l) and therefore the result follows from the 



fact that fj, is a probability measure. 



□ 



Remark 4.6 (An abstract result). It is clear that the proof of Theorem 4.5 is 
general and there is a similar statement in a purely abstract setting. The only 
two properties which we have used are (i) that the energy of is a weakly 
lower semi-continuous function of the first two density matrices; and (ii) that 
ejj(A) > e^(l) for all < A < 1. The first means that particles escaping to 
infinity always carry a non-negative energy. On the other hand, the second 
property means that, in Hartree theory, we can send some of the particles to 
infinity without any cost. The two assumptions are complementary and both 
are necessary here. By adding a constant to T, we can always ensure either (i) 
or (ii), but not both at the same time. 



Before proving Proposition |4.4| under the general assumptions of Lemma 4.3 
let us remark that this proof is obvious if w > 0. In this case (4.4) follows 



immediately from Fatou's lemma for operators and from the fact that the es- 
sential spectrum of T starts at 0, that is, 1(T < 0) is a compact operator. An 
example of such purely repulsive systems is a bosonic atom. 

Example 4.7 (Bosonic atoms). 

After scaling, an atom with a classical nucleus at the origin in M 3 and N 
"bosonic electrons" is described by the Hamiltonian 



JV 



1 



t\X 4 



+ 



1 



N 



E 



1 



1 1 — ' \Xh 

l<k<£<N 1 K 



(4.8) 



acting on Sj N = ® s L 2 (R 3 ), where t = (N — 1)/Z is the ratio between the 
number of ele ctron s and t h e nuc lear charge. This Hamiltonian satisfies all our 
assumptions ( | 1 . 16| )-( | 1 .1 7] )— d 1. 18| ) with d = 3, s = 1, V(x) = -(ijxj)- 1 and 
w(x) = \x\ 



Since w > 0, the equivalent assumptions of Lemma 



satisfied and (4.4) follows immediately. 



4.3 



are also 



For any fixed t > 0, the Hartree energy e^(A) is decreasing on [0, min(l, t c /t] 
and constant on [t c /t, oo), where t c ~ 1.21 018]. Furthermore, e^(A) admits 
a unique minimizer kq if and only if A < t c /t. We deduce that the set M. v 



appearing in the statement of Theorem 4.5 is 



M 



v 



{e i9 u } 



6»e[0,27r) 



where 



1 if t < t c , 
t c /t if t > t c , 
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■,V(-t\ ;f + s + „„J f^,. „V I 



where uq is the unique minimizer for e£(l) if t < t c and for e£(* c /i) if t > t c . 
Now we deduce from Theorem 14.51 that 



(fc) 
Tn 



: )(4 



strongly in & 1 (Sj k ), if t < t c , 
k - . i«i?*)(u? fc | weakly-* in 6 x (ij fe ), if i > t c , 

for all > 1. To our knowledge, this is the first result about approximate 
minimizers for bosonic atoms beyond the critical value t c . 

The occurrence of complete Bose-Einstein condensation was shown by Ben- 
guria and Lieb in [9] for t < t c (see also [62j |6]). Their proof used the 
particular form of the Coulomb interaction through the Lieb-Oxford inequality. 
Recently, Kiessling also considered bosonic atoms in |33] for which he used the 
classical de Finetti Theorem. 

We end this section by providing the 



Proof of Proposition ^ -4\ If TrsI-^Tjv''] ^ s n °t bounded then there is nothing to 
prove, so we may assume that Tr^j^yjy ] < C. The proof is divided into two 
steps. 

Step 1. Splitting of the energy. The first step is a classical result used 
many times in the literature, and which does not require that w has no bound 
state at infinity. 

Lemma 4.8 (Splitting of the energy). 

Assume that (1.16), ( |1.17| ) and (1.18) are satisfied. Consider a smooth partition 
of unity x 2 +?/=l where x{ x ) = 1 tf\ x \ ^ 1 an d x( x ) = tf x — 2, and d efine 



Xr{ x ) = x( x / 'R) an d VR = r](x/R). Let 7^ and jff be as in Proposition 4-4 
We have 



liminf ( Tr«[T7 



+ -Tr 



(2)- 



> liminf liminf ( Try-, [TxRjff Xr] + 7> Tr fl 2 ^Xr^n' Xr 2 

R— >oo TV— >oo \ Z 



92 (2)®2i 



+ Tr Si [Kr, R7 ^r lR ] + ^ Tr^wr,^ " V 

.,(1)1 



n ®2 (2) <g,2- 



(4.9) 



Proof of Lemma 4-8- Again we may assume that Tt^[Kj n ] < C. To split the 
kinetic energy, we use the estimate 

lim \\K - xrKxr ~ m K m\\L^L 2 = ( 4 - 10 ) 

R— >OD 

which follows for instance from the fractional IMS formula which can be found 
in |35| Lemma 7]. Moreover, using that K 1 / 2 "/^ K 1 / 2 is bounded in the trace- 
class, we have 



Tr(Vi] R i$r} R ) 



< C 



(1 - A)- s / 2 m V m (l - A) 



-s/2 



(4.11) 



Our assumption (1.17) on V implies that the norm on the right side goes to 
when R — > 00 and we obtain 

lim liminf Tr(Vr) R r y^rjji) = 0. 
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This means that the particles at infinity will never see the potential V. 

To split the two-body term, we insert 1 = (Xr( x ) + r /^( x ))(Xij(y) + Vniv)) 
and expand. We only have to control the cross term, namely we have to prove 
that (possibly for a subsequence) 



lim lim 

R— too N—too 



Xr(x)\w(x - y)\rj R (y)p i ^\x,y)dx dy 







(4.12) 



where p%\x,y) = i$(x, y; x, y). To this end we write T] R (y) = r] R (y)-r]j R (y) + 
Vht(y) and remark that Xr(x)\w(x - y)\r]j R (y) < l\\ x - y\>R\ w(x - y) which can 
be easily controlled using the same argument as for ( 4.11[ ). So it remains to 
treat the term with Xii( x )[Vit{y) ~ ^IrCs/)]) f° r which we claim that 



lim lim 

R— >oo k—^oo 



w(x - y)\x 2 R(x)[r]R(y) - ■nlR{y)]pkS x ^ v) dxd v 

< lim lim ( // \w(x — y)\t(R < \y\ < 8R)p^ (x, y) dxdy 

R— >oo k— >oo \J J k 







for an appropriate subsequence (iVfc). This is an adaption of Lions' concentration- 
compactness argument. Let us introduce the concentration functions 

Q N {R) := JJ \w(x-y)\l(\y\>R)p%\x,y)dxdy. 

For every N, the function R i— > Qn(R) is non-increasing on [0,oo). Moreover, 

< Q N (R) < TTtf[\w\^ ] ] < C by ( [TTT8| ). Therefore, by Helly's selection 
principle, there exists a subsequence and a decreasing function Q : [0, oo) — > 
[0, Co] such that QN k (R) — > Q(R) for all R £ [0, oo). Since liniR^oo Q(R) exists, 
we conclude that 



lim lim (Q N JR) - Q Nk (8R)) = lim (Q( R ) - Q&R)) = 0- 

R— >oo k— »oo R— >oo 

This is the desired convergence. 



□ 



Step 2. Passing to the limit. Now we look at the right side of (4.9). By 
the local compactness we have 



lim xflPxf = xfl^xf 

iv — >oo 



(4.13) 



strongly in the trace class for all k > 1 and any fixed R. Using the fact that for 
any one-body operator < A < 1 



AyffA - Tr 2 [A^>A^ 2 } >0onS), 



2„,(2) 



(4.14) 



(where Tr2 is the partial trace in the second variable) and Fatou's lemma, we 
get 



lim inf Try- 



(T + C ) (XR^XR-T^xf^xf] 



> Tr, 



(T + C^ixRl^XR-T^ixTl^xt 



2^(2). 



(4.15) 
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for a large enough Co chosen such that T + Co > 0. By Fatou's lemma again, 



lim inf - Tr 

iV-s-oo 2 



ft 2 



(T®1 + 1(8>T + 2Co + ^xf^xf 2 ] 
1 



> -Tr f 



(T®l + l»T + 2Co + w)x| 2 7 (2) xl"] • (4.16) 



Here we have to ch oose C p eve n larg er to make sure that (T®1-\-1®T+2Cq+w > 
on fj 2 . Putting (4.15) and (4.16) together and using the strong convergence 



(4.13), we find that 



liminf (^[TxrtWxr] + \ ^^xfl^xf 



> T^[T X r1 W Xr] + -Tr„4wxf7 i2) Xp 



(4.17) 



By taking now the limit R — > oo we arrive at 

lim inf lim inf ( Tt^TxrJ^Xr} + \ Tr fl2 [wx® 2 ^ Xr 2 ] J 

yTr^T^j + ^Tr^w^]. (4.18) 

Up to now the argument is general and it holds without the assumption that 
w has no bound state. We will only use this fact to estimate the second term 



on the right side of (4.9), that is, to prove that the particles far away have a 



non negative energy. Indeed, using (4.14) and the assumption that H°(2) > 0, 
we obtain 



Tr fi [Kruaftvd + \ ^ \^fl%\f\ > \ ^ ^vf^vf -0.(1.1 9) 



The weak lower semi continuity (4.4) now follows from (4.9), (4.18) and (4.18). 



□ 

4.3. Translation invariant case. In this subsection, we ignore the potential 
V and consider the ground state energy E°(N) of the fully translat ion-invariant 
Hamiltonian 

N , N 



A,; 



m 



+ 



N 



i<j 



We work under the same assumptions (1.16), (1.17) and (1.18) as before, but 



do not assume any of the equivalent statements of Lemma |4.3| The particles 
escaping to infinity can thus have a negative energy, and the simple proof of 
Theorem |4.5| in the previous section does not apply. Furthermore, the situation 
is complicated by the fact that we do not expect the local convergence of all the 
sequences of approximate ground states. By using the translation invariance 
and the kinetic energy of the center of mass, one can construct sequences for 
which 7^ (• — xn) ^* for any translation (xjv) C M. d . This is called vanishing 
in Lions' terminology }47j . Our proof of the validity of Hartree's theory will be 
based on some ideas of Lieb, Thirring and Yau [451 14"6] and on the geometric 
techniques of [36 j. Our main result is the following. 
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Theorem 4.9 (Translation-invariant systems). 

Under Assumptions (1.16), (1.17) and (1.18) onw, we have 



We recall that the ground state energy in Hartree theory is 

1 



<&(!):= inf Uu,Ku) + 



\\u\r=i 

where K := (m 2 — A) 6 



\u(x)\ w(x — y)\u(y)\ dxdy 



m 



2.s 



Example 4.10 (Boson stars). 

A pseudo-relativistic model of a star with TV gravitating bosons may be de- 
scribed by the Hamiltonian in dimension d = 3 with m = l,s = 1/2 and 
w(x) = — for some < k < 2/ir. The validity of Hartree's theory in 
this case was proved by Lieb and Yau |46j . Their proof is based on a clever 
replacement of the two-body potential by a one-body potential, which is also 
crucial in our proof. 



The outline of our proof is as follows. Following Lieb and Yau |46j . we use 
some part of the two-body potential to create a negative one-body potential 
which breaks the translation invariance of the system. Since e^(l) < by 
Lemma 4.2 it suffices to consider the case when limjv->oo E(N)/N < 0. In 



this case, we may rule out both the vanishing and dichotomy for the modified 
model, and hence the geometric method of |36] applies. The final result then 
follows from an approximation argument. 



Proof of Theorem J±.9. First, we remove the center of mass to create a negative 
one-body potential. To be precise, for any TV-body wave function ^ 6 Sj N and 
e > 0, we can write, similarly as in |45| 146] . 



N — 1 



N 



N-l 



N-l 



*' ( X) ( Ki ~ £W -( X i ~ X n)) + N \ 2 ^ W ^ Xi ~ Xj ^ 
\ i=l i<j 



(4.20) 



where w-{x) := max{0, — w(x)} and w e (x) := w(x) + 2eu>_(x). In (4.20), the 
Hamiltonian in the parenthesis depends on xn but, by translation-invariance 
of the other terms, the bottom of its spectrum is actually independent of xn- 
Therefore, 

E°(N) > E £ (N-l) 



N ~ N-l 

where E £ {N — 1) is the ground state energy of the modified Hamiltonian 
N-l ^ N-l 

H £ ,N-1 = ^2 { Ki ~ £w -( x i)) + N _ 2 ^2 W£ ( Xi ~ x ^ 



1=1 
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4.1 



inij N 1 . Since E°(N)/N and E £ (N)/N are increasing sequences by Lemma 
the following limits exist: 

EJN) , E°(N) n , . 

a £ := lim £ ^ T ; < lim — =: a < e&(l) < 0. 

We may now assume that a < 0; otherwise a = e^(l) = and the proof is 
finished. Note that 

lim inf a(K — ewJ) = inf a(K) = 0, 



where we have used the lower bound a~K + w > — C by (1.18). Therefore, 
inf o"(AT — ew_) > a for e > small enough. In this case we shall show that a £ 
is exactly equal to the modified Hartree energy 



e H , e (l) := inf < (u, (K - ew-)u) + 



1 



\u(x)\ 2 w e (x - y)\u(y)\ 2 dxdy 



Since a £ — > a and eH, £ (l) — >■ as e — >• 0, we then conclude that a = e^(l) 

as desired. 

In order to prove that a £ = eH, e (l), we consider a sequence of wave 
functions such that (^f N , H £jN ^ N ) = E e (N) + o(N). Then 



lim 

AM-oo 

„( fc ) 



N 



lim (t^[(K - ew-)y$>] + \ IVKtJ?]) 

A<— »oo \ Z / 



where 7]^ are the fc-particle density matrices of ^at. Up to extraction of a 

subsequence we may assume that 7)^ ^* 7^ for all k > 1. The equality 
a £ = eH, e will follow immediately if we can show that 



Tr^ 1 )] = 1. 



(4.21) 



Indeed, by Corollary 2.4, (4.21) implies that 7^ — > strongly for all k > 1. 
Then, arguing like in (3.10) for the proof of Theorem 3.1 in the confined case, 
one sees that, in this case, 



lim inf (tt^K - ew.)^] + \ Tr fi2 [w £ ^]) 

A<— >oo \ Z I 



The result then follows from the (strong) quantum de Finetti Theorem 2.1 



So it only remains to prove our claim (4.21 ). To this end we use the geometric 
localization method. We pick a smooth partition Xr+Vr = 1 like in Lemma 4.8 
and get 

a £ > lim inf lim inf { Tr^ [(K - ew^xrJ^Xr) + \ Tr fi 2 KxIS^Xr 2 ] 

R — >co A— >oo [ Z ... 

+ Tr fl [Kruaftw] + \ Tr# [w^f^rjf] } . (4.22) 
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We consider the xr~ an d ry^-localized states and CfL of defined 



in Q) and fl277| ). We have 

k 



Tr,[(/v - - ) X V ' x «] + i Kxf^xf] 

k 



and 



i=l i<j 



^N,k 



(4.23) 



1 N 

k=l 



k | k \ 

H ^ + ]\7^T ^ W ^ Xi - X ^ G N,k 



t=l i<i 

We have to estimate these two terms from below. 



(4.24) 



Estimate on (4.23). We apply the variational inequality 

A + tB = (1 - t)A + t(A + B) > (1 - t) inf <r(A) + i inf <r(A + 5) (4.25) 

for A = ELi( K - ew -)^ A + B = H £tk and t = (k - l)/(N - 1). Note that 
in this case inf cr(A) > inf a (A + -B) since 

inf a(K — ew-) > a > a £ > k^ 1 inf a(H £ ^). 



Thus from (4.23) it follows that 



rr UTS \ (!) 1 , 1 TV r 02 (2) ®2i ^ ^ ^ G iV,fc -^eCO 

Tr^K - ew_)xR7]v Xfl] + ^ ^VKx^ 7jv ] > 2_> jy F~ ' 



fe=l 



Since 



^[XbtSJ ] ,7"^ Tr [x|7 (1) ] and lim 



fc=0 

we conclude that 



k— >oo k 



lirninf (Tr^[(K - ew^xr^xr] + \ ^{w^f^xf]} 

>a £ Tr[ X | 7 {1) ]- (4.26) 



Estimate on ( |4.24[ ). Using if > 0, w £ = w + 2ew_ > (1 - 2e)u> and 2?°(&) < 
a/c < 0, we find that 



1=1 



(l-2e)(fc-l) Q 
iV^l ^ 



> (1 2e){k l) E\k) >E°(k)-2eak 
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for every k > 1. Thus from (4.24) it follows that 



rr \ts (!) 1 , 1 t r ®2 (2) 02i ^ k Tr G N,k ( E°(k) \ 
Tr^KrjR^'m] + - Tr fi2 [w e ri% 7jv >Vr ] > ^ ^ ( — r 2eo 1 . 

fc=i ^ ' 



Since 



N 



E 



Tr 



N,k 



TtIVrIn) 1 " Tr[x|7 (1) ] and lim 



fc=0 

we deduce that 



liminf (TtssFviaWw] + l^fA^rifl^vf}) > (l-2 £ )a(l-Tr[ X | 7 (1) ])- 

(4.27) 

Substituting now ( |4.26| ) and ( |4.27[ ) into ( |4.22[ ), we find that 

a £ > liminf (o e TV[ X |7 {1) ] + (1 - 2e)a(l - Tr[ X | 7 (1) ]) 
= a £ Tr[ 7 «] + (1 - 2e)a(l - Tr[ 7 «]) 



Since a e < a < 0, we conclude that Tr^ 1 )] = 1 as stated in (4.21) and the 
proof is complete. □ 



4.4. General case: Proof of Theorem 1.1 We are now able to prove our 
main result, Theorem |1.1[ which was stated in the introduction. Our strategy 
is to split the energy into two parts corresponding to the particles staying in a 
neighborhood of and those escaping to infinity. We use the weak de Finetti 



Theorem 2.2 for the local part. The particles far from the origin form a fully 
translation-invariant system for which we have already shown that Hartree's 
theory is valid. The conclusion then follows from the binding inequality in 
Hartree's theory. 

Proof. Let V^jv be a sequence of wave functions such that n , jy) = 

E V (N) + o(N) and denote by 7 jy the &;-body density matrix of *&n- Up to 

a subsequence we may assume that 7^ — 7^) weakly-* in the trace class for 
all k > 1. Let us denote by \i the probabilit y me asure on BS) associated with 
1 as in the weak de Finetti Theorem 



2.2 



First we proceed similarly as before. Let Xr + f] R = 1 be a smooth partition 
of unity as in Lemma 



4.8 



G^y in the Fock space By (4.9), we have 



and consider the associated localized states and 



lim 

TV-s-oo 



E V {N) 
N 



lim 

N->oo 



(2)- 
N - 



> liminf liminf \ Tt^Txr^ Xr] + \ Tr^xfl^xf] 
+ Tr fi [KmaffvR] + \ [vmfl%\T] } • (4-28) 
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By the strong local compactness (as in (4.17)) and the weak quantum de Finetti 
Theorem 12.21 we infer 



liminf iTr^TxRlffxR) + \ ^ [wxfl^xf]] 



£%{XRu)dn{u). (4.29) 



For the second term of the right side of (4.28), by using the geometric localiza- 
tion method we can show that 



.X 1 ) 



1 



„®2 „,( 2 ) <X>2- 



liminf Tr[Tr] R f N ) r] R ] + -Tr^fwyg 1 N Vr 

iV— >-00 1 



> 



e H (! " \\XRu\\ )dfi(u). (4.30) 



BS) 



Before proving (4.30), we explain how to conclude the proof of the theorem. 
By substituting (4.29) and (4.30) into (4.28) and using Fatou's lemma we find 
that 



lim 

AT-s>oo 



E V (N) 
N 



> lim inf 

R— >oo 



£K(XRu) + e^(l 



\XRU\ 



> 



> 



Bfj 



BSj 



BSj 



lim inf 

R-^-oo 



£r(xru) + e^(l - \\xru\\ 2 ) d/j,{u) 



£g(u) + 4(l 



dfj,(u) 



e^CII^H 2 ) + egj(l — H^ll 2 )] d(i(u) > e H (l). 



Here we have used the continuity of A h-» e^(A), which will be proved be- 
low, an d th e large inequality e^(l) < e^(A) + e^(l — A) which is taken from 
Given the upper bound E V {N)/N < e£(l), we conclude that 



Lemma 



4.2 



limjv-j.oo E(N)/N = e^(l) and that \i has its support in 

M V = {u G BSi : Sl{u) = eUhf) = e£(l) - 4(1 - ||n|| 2 )} . 
Moreover, when e£(l) < eg(A) + e^(l - A) for all < A < 1, we can deduce 



stronger statements as in the proof of Theorem 3.1 



Now we prove (4.30). Let us consider the ^-localized state G V N of which 
is such that 



1 



92^(2). 



1 N 

jvE^V 



k=l 



\i=i i<j j 



k=l 



N TrG^ 



N 



Wi 



i=l 
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On the other hand, by the fundamental relation TrG^ fc = Tr N _ k (men- 
tioned before in (2.8)) and Theorem 2.6 we have 



N 



lim ^TrG^e° H 



A:=0 



N 



N 



lim ^TrG^_ fc e° 



fc=0 



N 



N / h\ 

lim YTrG*. el (1- — J 
k=0 v 7 



e^ 1 - \\XRu\\ 2 )d^(u)- 



Therefore, in order to prove (4.30) it suffices to show that 

JV 



lim YTrG v Nk 

fc=0 

where &i(A) = and 



k 



k-1 



N \N - 1 



o, 



(4.31) 



1 ( k X k \ 

b k {\) := - inf o& K t + j—^ ^ Wij 

\i=i i<j J 



when k > 2. 



Note that we have to deal here with a Hamiltonian of the same form as but 
with a factor A € [0, 1] in front of the interaction. By Theorem 4.9, we know 
that Hartree's theory is correct for such Hamiltonians, that is, we know that 



lim Xb k {X) = A inf ( (u,Ku) + 



k— >oo 



Nr=i 



A 



iMr=A 



inf (u, Ku) + 



1 



w{x - y)\u{x)\ 2 \u(y)\ 2 
w(x-y)\u(x)\ 2 \u(y)\ 2 )=e° }1 (X). 



So we are almost done. In order to justify (4.31), let us show that the func- 
tions {frfcjfcLi are equicontinuous on [0, 1]. By adapting the variational estimate 
(4.25), we obtain 

6fc(A) > bfc(A') for all < A < A' < 1. 
On the other hand, if we denote 5 := (A 7 - A) (a -1 - A)" 1 , then 



11 A' 

\i=l i<j 



1-5 



+ 



ka 



i=l i<j / 

\ »=1 i<j J 



> (l-S)b k (X) 



C5 



a 



Thus 



< b k (X) - b k (X') < S(b k (X) + Ca- 1 ) < C\X' - X\ 
for a constant C > independent of A, A' and k. 
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The equicontinuity of {bk}^ =1 and the pointwise convergence limfc_ ) . 00 A6fc(A) 



e^(A) yield the uniform convergence 



lim sup 

M-too N>k>M 

Consequently, we find that 



lim sup 

N^oo k=l,2,...,N 



-b ( — 

N h \ N 



-b ( — 
N h \N 



N 

k 
N 



0. 



and (4.31) follows. The proof is complete. 



□ 



5. Further extensions 



We conclude this paper by mentioning four interesting cases that may also 
be dealt with using our method. 

Remark 5.1 (Bosons in a magnetic field). Our results in Theorem |1.1| are still 
valid when the fractional Laplacian (m 2 — A) s is replaced by its magnetic version 
(m 2 + \X7 + iA(x)\ 2 ) s , where A : ]R — > M rf is a Borel measurable ve ctor p otential. 
For simplicity we assume that |^4| 2s = fs + fe with fj being as in (1.17). In this 
case, the IMS-type estimate 

\K - XrKxr ~ r]RKr] R \\ L 2^ L 2 = (5.1) 



lim 



in (4.10) still holds true with K replaced by (m 2 + |V + iA(x)\ 2 ) s . The proof 
of (5.1) follows the same argument as [35., Lemma 7]. Moreover, by using 
the Cauchy-Schwarz inequality and the operator monotonicity of t \— > t s when 
< s < 1, we can show that 



;i-a) 



-*/2. 



Vr[ { m2 + |V + iA(x 



\2\s 



(m 2 -A) s W(l-A) 



-s/2 



0, 



lim 

which is a substitution for (4.11). Therefore, our approach applies exactly as 
in the non-magnetic case. 

Remark 5.2 (Bosons hoping on a lattice). In this paper we mainly considered 
continuous systems for simplicity. Our method applies as well to bosons living 
on a lattice £ C M rf (a discrete subgroup of W 1 with compact fundamental 
domain) and with a kinetic energy described by the discrete Laplacian. In this 
case we simply assume that the potentials V and w are in £°°(>C) and tend to 
zero at infinity (they are then compact operators on Sj = £ 2 (C)). As there is 
always local compactness on the lattice and as the discrete Laplacian satisfies an 



IMS localization formula similar to (5.1), our method applies mutatis mutandis 



and Theorem 1 1 . 1 1 holds in this well, without any change. 

Remark 5.3 (The absolute ground state). Our method may also be used to 
investigate the absolute ground state energy of a quantum mechanical Hamil- 
tonian, that is the infimum of the spectrum with no symmetry restriction on the 
admissible states. The absolute ground state energy coincides with the bosonic 
ground state energy in the situation covered by Theorem |1.1| by a well-known 
method |42l Section 3.24], but it need not be the case in general, for example in 



the presence of magnetic fields or rotation (see Example |3.4| and Remark 5.1). 
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Observe that any Hamiltonian Hn of the form (1.1) satisfies 



Uct n HnU* n = Hn 

for all permutations ctn of N variables, where U aN is the unitary operator 
permuting the order of variables according to cjm- Therefore, when analyzing 
the absolute ground state energy of fljy we may consider only the mixed sym- 
metric states, which are positive trace-class operators acting on $j® N and 
satisfying 

^7 W f/;=7 (W) (5.2) 
for all permutations ctjv- In this language the assumption on the Bose-Einstein 
symmetry corresponds to the stronger condition 

^7 W =7 W ^=7 W (5.3) 

for all permutations on- 

If we are given an infinite sequence {7^} fc _ of fc-particle positive trace- 



class operators that satisfy the symmetry assumption (5.2) and the consistency 



assumption (2.1), then the quantum de Finetti theorem [31J states that there 
exists a Borel probability measure \i on the unit sphere S6 (Sj) of trace class 
operators, invariant under the group action of S , such that 

7 (fc) =/ 7®*rfM(7) (5-4) 

for all k > 0. It is then not difficult to adapt our approach to deduce a weak 
version as in Section [2] and use it to prove results about the absolute ground 
state that parallel those we presented for the bosonic ground state. 

In this case, one obtains as limit object a Hartree theory for mixed one-body 
states, that is the ground state energy per particle is given in the limit N — > 00 
by the minimization of the functional 

£ H ( 7 ) == Tr^[r 7 ] + l - TV^a [u> 7 ® 2 ] (5.5) 
over all positive trace-class operators 7 on fj with Tr 7 = 1. 



The general question of when the minimization of (5.5) reduces to that of 



( 1.2 ) (that is when absolute minimizers are asymptotically bosonic in the mean- 
field limit) seems to be mostly open. In [55] the absolute ground state energy 
of a rotating trapped Bose gas with repulsive interactions has been considered 



in the Gross-Pitaevskii limit, and a functional similar to (5.5) has been derived. 



In this particular case, sufficient and necessary conditions are also given for 



the minimization of (5.5) to reduce to that of (1.2). These conditions are inti- 
mately linked to the question of symmetry breaking and nucleation of vortices 
in rotating Bose gases. 

Appendix A. Alternative proof of the weak de Finetti theorem 



Here we explain how to derive Theorem |2.2| from the strong de Finetti theo- 
rem, following ideas of Hudson and Moody [31] . This proof is more direct but 
not as constructive as the one we gave in Section [2j which also allowed us to 



relate the de Finetti measure \x to geometric localization in Theorem 2.6 
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Let us consider a sequence of normal states G 6 1 (Sj N ) and the weak-* 
limits y( k \ k > 1, of their reduced density matrices, as in the statement of 



Theorem 2.2 Following [52{ [70], it is useful to think of Tjv as a state wtv on 
B(f)®°°) which is the C* inductive limit of the sequence B{Sj® N ) of all bounded 
operators on the tensor product Sj® without symmetry. This means that 

wtv(&i ® . . . ® 6m) := Tr fl ®iv (T^ b\® ...® b N ) 

for all M > N and where T^v is extended to outside of S) N C fi® N . Here we 
have made the abuse of notation to identify the operator b\ ® . . . ® bu with 
b\® . . .® bu ® 1 &> — Note that, using the bosonic symmetry of 

ujNih ® ■■■® bu) = ^N(b a (i) &) . . . &) 6 ct (at) 8) 6at+i <8> • • • <8> b M ) (A.l) 

for any M > N and any permutation cr of the first N variables. 

By the Banach-Alaoglu Theorem, the sequence ujn admits a weak-* cluster 
point uj, a state on B($j®°°). So ujn converges to uj along a subnet, which means 
that for any n and any bi,...,b n G B(S)) 

Uh(a) {h ® ■ ■ ■ ® b n ) ->■ uj(bi <g) . . . ® b n ) (A. 2) 

where h : A i— > N is a monotone cofinal function from some directed set A to 
the integers. It is of course important to be able to test against the identity 



operator in (A. 2), to ensure that a; is a state. 

The state uj determines a hierarchy of n-particle states uj( n ' on B(!f)® n ) with 
the consistency relations 

J m \bx® . . .®b m ) = uj(b 1 ®...®b m ) = J n \b 1 ®...®b m ®...®l® n - m ) (A.3) 

for all n > m, and all bi,...,b m 6 -6 (•£>)• 



From (A.l) and (A. 2) we deduce that the cluster point uj is symmetric, that 
is, uj^(b a m ® ■ ■ • <g) &o-( ra )) = u/ n )(&i <g> • • • <8> b n ) for all bounded operators 
bi,...,b n G £>(f)) and every permutation u. Now, from the strong de Finetti 
theorem (for abstract states on an algebra) of |65| 131] . there exists a Borel 
probability measure /i on the set of states S($)) on S) such that, for any n > 0, 

= / dp^d*". (A.4) 

This is a consequence of the fact, proved first in |65j . that the tensor powers 
uj®°° are the extreme points of the convex set of symmetric states on B(Sj°°). 
The link between uj^ and the weak-* limits 7^ of the density matrices of the 
sequence T^v is that 



Tr.iv i^K) = lim Tr.iv K) = lim uj n (K) = uj(K) = uj (n) (K) 



(A.5) 

for every symmetric compact operator K on Sj N . 

We recall that any state uj on B(&) can be restricted to the algebra K,{3) of 
compact operators (M is here any fixed separable Hilbert space). The associated 
non-negative linear form uj nor is called the normal part of uj and it necessarily 
arises from a trace-class density operator 7^ £ S 1 (^), since 6 1 (^) is the dual 
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of tC(&). Indeed, the map u h- > 7^ is continuous from the set of states into the 
trace-class and we have 

Tr(7 w ) = sup Tr(7 w K) = sup u(K) < sup u(B) = oj(1) = 1. 

0<K<1 0<K<1 0<B<1 



By (A.5), we deduce that 7^) must be the density operator associated with the 



normal part of the state u/ n ) . 

The main point of the proof is the remark that the normal part of a tensor 
product is the tensor product of the normal part, 

(c^) nor = ( Wnor r, (A.6) 
which follows from the fact that !C(Si® N ) = JC(fi)® N , From this we deduce that 

7 (n) = ( d^)( lul )® n . (A.7) 

Since, as we have said above, the map u 1— )■ 7^ is continuous, we can consider 
the push- forward Borel probability measure jl on the unit ball B& 1 (S)) of the 
trace-class, which is such that jt(A) = /i({u> : 7^ 6 A}) for all Borel sets 
A C Be 1 (S)). Hence 

7 {n) =/ dmi® n - 

JB6 1 (^i) 

To conclude the proof of the weak de Finetti theorem, there only remains to 
show that jt is supported on the set of pure states: jt ({\u) (u\ , u € B$S}) = 1. 
This follows exactly |31| Section 4] and this is where we need the important fact 
that our state has the Bose-Einstein symmetry, that is, S n ^ n ' = 7^ where S n 
is the symmetrization operator. Taking the trace against S n we find 



Tr 7 W = Tr S n 7 W = f dp.{rt) Tr (S„ 7 
JBe 1 (h) 



< / dfii-y) Tr 7 ® n = Tr 7 ( n ) 
JBG 1 ^) 

where we have used that S n < 1. From this we deduce that Tr(S n ~/® n ) = 
Tr 7 ® n for all n > 1 and /2-almost all 7. This is equivalent to 7 = |it)(tt| 
by [3H Proposition 3]. Therefore jl is supported on rank-one density operators. 
Associated with this measure, there is a unique ^-invariant Borel probability 
measure // on the ball Bfy such that 

7 H= f d^{u)\u® n ){u® n \, 
and /u' is the sought-after de Finetti measure of Theorem 2.2 □ 
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